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A bstract 

The CCSD(T) method has been used to compute a highly accurate quartic force 
field and fundamental frequencies for all 14 N and 15 N isotopomers of the high energy 
density material T d N 4 . The computed fundamental frequencies show beyond doubt tha 
the bands observed in a matrix isolation experiment by Radziszewski and coworkers are 
not due to different isotopomers of T d N 4 . The most sophisticated thermochemical 
calculations to date yield a N 4 2N 2 heat of reaction of 1 82.22±0.5 kcal/mol at OK 
(180.64±0.5 at 298K). It is hoped that the data reported herein will aid in the ultimate 

detection of T d N 4 . 

Introduction . 

Beginning in the early 1 990s, theoretical quantum mechanical calculations 

identified tetrahedral (T d ) N 4 as a potentially useful novel rocket fuel for a number of 
reasons [1-11]. First, a considerable amount of energy is released upon dissociation into 
two N 2 molecules ( 1 83 kcal/mol [1 ]; see also Ref. 2 in [12]). Second, there is a large 
barrier to dissociation into two N 2 molecules, and the curve crossing to the lowest triplet 
state was 13 kcal/mol above the T rf N 4 minimum making it possible to handle T d N 4 . 

Third the main dissociation product, N 2 , is environmentally benign. Further studies 
examined the spin-orbit crossing probability as well as the barrier to dissociation, and 
concluded that indeed T d N 4 held promise as a novel rocket fuel should it be possible to 

synthesize the compound. . . 

Recently, there has been renewed interest in the high energy density material 1 d 

N, [12-15], Aside from the study by Radziszewski and coworkers, which will be 
discussed below, the excited electronic states of r rf N 4 were investigated [12,14], In one 
of these studies a novel approach to the synthesis of T d N 4 was proposed [14], whic is 
based on combining two N 2 molecules in bound quintet excited electronic states. i e 
this new approach appears promising and is the first proposed method which targets the 
electronic structure of the reactants to be similar to the desired end result, T d N 4 , it has 
yet to be attempted experimentally. In addition, Cacace et al. [15] have very' recently 
reported the synthesis of a tetranitrogen compound, although this is not T d N 4 . Further, 
this compound is not a high energy density material, and since it has a chain- like structure 
with two short N-N bonds, it is difficult to see how it will be useful in making T d 4 . 
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In the summer of 2000 it appeared that at last some progress had been made when 
a study was published by Radziszewski and coworkers in which they tentatively claimed 
to have made 7^N 4 [13]. The claim was based on observation of the one infrared active 
band of Tj N 4 and its full 15 N isotopomer. However, the isotopic shift of the band did 
not agree with theoretical calculations with the experimental value being 36.7 cm while 
theory placed this shift at 31.6 cm' 1 [13,14], Radziszewski and coworkers were also 
concerned about their inability to observe IR bands from any of the mixed isotopomers. 
However, as we have recently explained [14], this is not surprising since observation of 
the 14 N 4 and 15 N 4 isotopomers was at their instrumentation limits, and the IR bands tor 
the mixed isotopomers would have significantly reduced intensities due to a splitting of 
the triply degenerate vibration. Thus the main concern from the work m Ref. [13] is the 
disgreement between theory and experiment for the 14 N 4 to 15 N 4 isotopic shift. Since 
matrix isolation was used in Ref. [13], it is possible to have bands shift due to mteracfton 
with the matrix, but it is highly unlikely that any shift would be different for N and N 
isotopomers, so this possible source of error for the isotopic shift can be eliminate . 
Another source of error, which is more of a possibility, is that the isotopic shift 
determined from theory is in error since only harmonic frequencies were used. We have 
addressed this issue in the present study by computing a highly accurate quartic orce 
field for 7^N 4 using the singles and doubles coupled-cluster method that inclu es a 
perturbational correction for triple excitations, CCSD(T) [16]. This method has been 
shown to yield highly accurate vibrational transition frequencies [17-20], and it should 
allow us to determine definitively whether the bands observed by Radziszewski and 
coworkers are indeed due to Tj N 4 . 

The theoretical approach is given in the next section while the results and 
discussion are presented in the following section. Conclusions are given in the final 

section. 


Theoretical Approach 

The CCSD(T) electronic structure method was used to obtain a quartic force field 
for 7b N 4 . In these calculations, Dunning’s [21] correlation consistent basis set that 
includes up through g - functions on nitrogen, cc-pVQZ, was used. In order to account or 
the effects of core-correlation, a CCSD(T) quadratic force field was obtamed with Woon 
and Dunning’s [22] cc-pCVQZ basis set. In all CCSD(T) calculations, only sphenca 
harmonic components of d,f and g-functions were used. As noted previously, the 
CCSDfT) method has been shown to yield highly accurate force fields and vibrational 
frequencies for a range of molecular systems [17-20], All CCSD(T) electronic structure 
calculations have been performed with the MOLPRO program package [23] . 

The equilibrium structure of the molecule was determined such that the residual 
internal coordinate gradient was less than lxlO* 6 atomic units. The quartic force field was 
then determined in symmetry internal coordinates. These coordinates are given by, 


S^aO = 1/V6 (r, + r 2 + r 3 + r 4 + r 5 + r 6 ), 
S 2a (e) = 1/Vl2 (2ri - r 2 - r 3 + 2r 4 - r 5 - r 6 ), 
S 2b (e) = l^ 4 0*2 - r 3 + r 5 - r 6)> 


( 1 ) 

( 2 ) 

(3) 


(4) 

(5) 

( 6 ) 


S 3x (t 2 ) = 1/V2 (r 2 - r 5 ), 

S 3 y (t 2 ) = 1 /V 2 (r 3 - r 6 ), 

S 3z (t 2 ) = 1/V2 (r x - r 3 ). 

The numbering of bonds is such that if r, connects one pair of atoms then r 4 connects the 
other pair and similarly for the {r 2 , r 5 } and {r 3 , r 6 } pairs. These are the same symmetry 
coordinates used previously for a quartic force field calculation for the tetrahedral 
molecule Be 4 [24], except that a typographical error for the S 2a coordinate has been 
corrected here. The quartic force field constants, given in Table 1, correspond to an 
unconstrained summation also as used previously in Ref. [24]. Only the symmetry 
unique force constants are listed in Table 1 — see Ref. [24] for symmetry relationships 
for the quadratic, cubic, and quartic force constants. Also presented in Table 1 are the 
quadratic force constants obtained at the CCSD(T)/cc-pCVQZ level of theory. The 
anharmonic spectroscopic analysis was performed with a modified version of the 
SPECTRO [25] program. 


Results and Discussion 

A. Anharmonic Constants and Accurate Vibrational Transition Energies 

The CCSD(T)/cc-pVQZ quartic force field given in Table 1 was used to obtain 
fundamental vibrational frequencies for the 14 N 4 , 15 N 4 , 14 N 3 15 N, l 4 N 15 N 3 , and 1 N 2 N 2 
isotopomers. The CCSD(T)/cc-pVQZ harmonic and fundamental frequencies are 
presented in Table 2. CCSD(T)/cc-pCVQZ harmonic frequencies for each isotopomer 
were obtained with the CCSD(T)/cc-pCVQZ quadratic force field given in Table 1 . Best 
Estimates for the gas-phase fundamental frequencies were then obtained by adding the 
CCSD(T)/cc-pVQZ anharmonic correction to the CCSD(T)/cc-pCVQZ harmonic 
frequency. These frequency values are also given for each isotopomer in Table 2. We 
note that the first two isotopomers have T d symmetry, the next two have C 3v symmetry, 
and the last one has C 2v symmetry. Elowever, we have ordered the modes for each 
isotopomer so that they are consistent with the parent l 4 N 4 isotopomer. 

Considering the two T d isotopomers first, it is apparent that core-correlation 
increases the frequencies by between 3.4 to 6.3 cm ' 1 for l 4 N 4 and by between 3.2 to 6.2 
cm ' 1 for 15 N 4 , which is somewhat less than what has been calculated for the N 2 molecule 
[26]. This is to be expected, however, since the N-N bonds in T d N 4 are single bonds 
whereas the N-N linkage in N 2 is a triple bond. Further the difference in the type of N-N 
bonds is also evident by some of the anharmonic corrections. For the N 4 isotopomer, 
anharmonic corrections for the two lowest energy modes are 16.9 and 18.4 cm-1, which is 
significantly less than that for the N 2 molecule, 28.6 cm-1 [27], However, the correction 
for the highest energy mode is 32.6 cm' 1 , which is probably a result of stretching (or 
compressing) all six bonds in T d N 4 simultaneously. Nonetheless, the t 2 mode is of most 

interest since it is the only one that is IR active [1]. 

Based on experience [17-20], we conservatively expect the Best Estimate 

fundamental frequencies to be accurate to within ±4cm ', with the Best Estimate 
probably being somewhat too high relative to the gas-phase value since the main 
remaining error is probably neglect of higher-order electron correlation effects beyond 



what is included at the CCSD(T) level of theory. Thus, taken in this light, our prediction 
for the t 2 mode for 14 N 4 of 935.3 cm' 1 would tend to support the experimental 
assignment, 936.7 cm' 1 , suggested in Ref. [13], keeping in mind that matrix effects can 
shift this somewhat as well. However, the residual error in our predicted isotopic shift 
for the !5 N 4 isotopomer, 30.9 cm' 1 , should be well under ±1 cm ‘ and this is not at all 
consistent with the experimental value of 36.7 cm ] . This leads to the interesting 
conclusion where we cannot definitively rule out an isotopomer of 7^N 4 for one of the 
bands observed in Ref. [13], but we can definitively say that both bands cannot be due to 
different isotopomers of T d N 4 . However, it seems unlikely that one of the bands is due 
to an isotopomer of T d N 4 and the other band due to a different as yet unassigned 
molecule. More likely is that the two bands belong to 14 N and 15 N isotopomers of the 
same molecule, but this molecule is most definitely not T d N 4 . 

Examination of the Best Estimate fundamental frequencies for the remaining 
isotopomers shows that of the bands originating from the t 2 mode in the 14 N 4 parent _ 
isotopomer, none of them are as low in frequency as the corresponding band in the j N 4 
isotopomer. In other words, the largest isotopic shift occurs between the N 4 and N 4 
isotopomers, so it is not possible that a mixed isotopomer accounts for the 
experimentally observed isotopic shift. From consideration of the data contained in Table 
2, there is little doubt that T d N 4 has yet to be synthesized and observed. 

For completeness, anharmonic constants according to the convention for spherical 
tops [28-30] are presented for 14 N 4 and ,5 N 4 in Table 3. It is hoped that the data 
presented in Tables 2 and 3 will be useful in the ultimate synthesis and detection of the 
high energy density material T d N 4 . 


B. Binding Energy Relative to Two N 2 Molecules 

In the decade since T d N 4 was first identified as a high energy density material, 
techniques for computing highly accurate thermochemistry have evolved significantly. In 
many cases it is possible to compute these quantities more accurately than can be 
obtained from experiment, with the exception of high resolution spectroscopy where this 
approach is applicable. The more sophisticated approaches, for example, take into 
account extrapolation to the one-particle basis set limit, core-valence correlation, and 
Darwin mass- velocity and spin-orbit relativistic corrections. In the present study we 
have utilized the W2 approach, which has been shown to yield atomization energies and 
heats of formation to better than ±0.5 kcal/mol for a set of 14 benchmark molecules [31], 
The interested reader is referred to Ref. [31] for details of the W2 approach. The data 
used in the W2 approach are compiled in Table 4 for the atomization reactions of N 2 and 
T d N 4 and for the dissociation of 7^N 4 into two N 2 molecules. CCSD(T)/cc-pVQZ 
reference geometries were used in these calculations. 

On examination of the data in Table 4, the first point to note is the excellent 
agreement between experiment and W2 theory for the N 2 molecule (being the reference 
compound, AH t;298 is 0.0 kcal/mol). This excellent agreement indicates that the W2 value 
for the heat of formation of 7^N 4 and for its dissociation into two N 2 molecules should be 
accurate to better than ±0.5 kcal/mol. Thus the important result from Table 4 is that the 


dissociation of 7^N 4 into two N 2 molecules is exothermic by 182.22±0.5 kcal/mol at OK 
or 180.64±0.5 kcal/mol at 298K. These values are in reasonable agreement with earlier 
calculations [1,12], but their uncertainty is smaller. Nonetheless, the important point is 
that they confirm beyond doubt that T rf N 4 is indeed a high energy density molecule. 


Conclusions 

A highly accurate quartic force field has been computed for the high energy 
density material T d N 4 . The CCSD(T) electron correlation technique has been used in 
conjunction with the cc-pVQZ and cc-pCVQZ basis sets. The CCSD(T) quartic force 
field has been used to compute accurate anharmonic corrections to the fundamental 
vibrations of T d N 4 including all 14 N and 15 N isotopomers. After correcting these for core- 
valence correlation, best estimate fundamental vibrational frequencies have been obtained. 
The isotopic shift between 14 N 4 and ,5 N 4 for the one IR active mode shows beyond doubt 
that the bands observed in a matrix isolation experiment [xx] cannot be due to different 
isotopomers of T d N 4 . It is hoped, however, that the improved fundamental vibrational 
frequency predictions for all of the 14 N and 15 N isotopomers of T d N 4 reported herein will 

be useful in the ultimate identification of T d N 4 . 

W2 theory has also been applied to the thermochemistry of 7^N 4 , allowing for 
the prediction of an improved dissociation energy with respect to two N 2 molecules. 
Specifically, W2 theory yields dissociation energies of 182.2210.5 kcal/mol and 
180.6410.5 kcal/mol at OK and 298K, respectively. These values show beyond doubt 
that T d N 4 is indeed a high energy density molecule. 
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table 1. LULHi;cq uiuuiiui 

li ii j unvi ^i.vv - 

CCSD(T)/cc-pVQZ 

CCSD(T)/cc-pC V QZ 1 

r„ 

1.4551 

1.4516 

F,. 

3.56221 t 

3.59671 

F-n 

4.42472 

4.46560 ^ 

F™ 

3.71350 

3.75208 

Fin 

-8.54584 



-11.12620 1 



-9.80063 



8.56000 1 



-9.36814 


1 -T 

-2.32956 


F 1 1 1 1 

18.24285 


r I 1 lota 

26.56394 


F 1 1 Ivtv 

22.73382 


r 1 1 o') o 

-9.19729 


r i 

27.34851 


r i 

5.00939 



68.73275 



31.34701 


F T k')h1'7'l'7 

-5.25058 



54.23396 


v 3.29649 

^ 3x3x3v3v 

T/X2 - T/ X 3 ^ T/ X 4 


Units tor quadratic, cudic, auu qumu^ * ? 7 

respectively. The equilibrium bond distance is in units of A. Only the unique symmetry 
internal coordinate force constants are given. See Ref. [xx] for symmetry relationships 
needed to complete the force field. 



Table 2. Theoretical harmonic, fundamental, and best estimate vibrational frequencies for 



in units of cm' 1 . The ordering of the modes is designed to be consistent with the parent 
isotopomer. 




r 


0.867 






Table 4. Thermochemical data (kcal/mol) for the noted reactions used in the W2 level of 


Liicuiy . 

2N — » N 2 

4N — » N 4 

2N2 — ^ N4 

SCFlimit- 

SCF/AV5Z 

0.04 

0.01 

-0.07 

CCSDlimit- 

CCSD/AV5Z 

1.64 

3.39 

0.12 

(T)limit-(T)/VQZ 

0.26 

0.77 

0.24 

SCF 

119.70 

8.50 

-230.91 

CCSD 

98.23 

238.41 

41.95 

(T) 

9.45 

28.02 

9.11 

CORE 

0.73 

0.31 

-1.14 

DMV 

-0.11 

-0.19 

0.04 

S-0 

0.00 

0.00 

0.00 

Equilibrium Total 

228.00 

275.05 

-180.95 

ZPVE 

3.36 

7.99 

1.27 

ZP Total 

224.64 

267.06 

-182.22 

AH 0 

0.42 

183.06 

182.22 

AFI298 

0.42 

181.48 

180.64 

t 1 1 


a See Ref. [xx] for a description of the W2 approach to thermochemical calculations. 




